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Abstract: Recently the non-extensive approach has been used in a variety of ways to describe dense nuclear matter. They differ 
in the methods of introducing the appropriate non-extensive single particle distributions inside a relativistic many-body system, 
in particular when one has to deal both with particles and antiparticles, as in the case of quark matter exemplified in the NJL 
approach. I present and discuss in detail the physical consequences of the methods used so far, which should be recognized before 
any physical conclusions can be reached from the results presented. 
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It is nowadays widely accepted that in many branches of physics experimen¬ 
tal data indicate the necessity of a departure from the standard extensive 
Boltzman-Gibs statistics, which is then replaced by a non-extensive statistics 
|T] . In high energy scattering one uses Tsallis statistics and Tsallis distributions 
which describe particle distributions over a large energy region with the help of 
only one single additional parameter, the non-extensivity q ^ 1. This replaces 
all additional parameters necessary when the extensive description is used; for 
q = 1 one recovers the usual BG statistics (cf. f2] and references therein). 
When considering dense Nuclear Matter (NM) some kind of mean held theory 
is used, either in terms of nuclear degrees of freedom (like the Walecka model 
(WM) [3]), or in terms of quark and anti-quark degrees of freedom (like in the 
Nambu-Jona-Lasino (NJL) model |4ll5ll6lj ). Both descriptions has been reformu¬ 
lated using a non-extensive approach (for the g-WM in [ I71ISMID] and for the 
g-NJL model in p0. All these models try to find the non-extensive traces in 
the nuclear Equation of State (EOS) using the non-extensive single particle dis¬ 
tributions. However, it is very difficult to compare their results in a conclusive 
way because they use different forms of single, non-extensive particle distribu- 

* Other approaches to non-extensive dense matter can be found in |12ll3ll4| . 
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tions. In this work, analyzing a Fermi Gas model we compare and discuss these 
distributions following the method presented in [IB], both for the momentum 
dependent fermion distributions n q (p) and for anti-fernrion distributions n q {p). 


In the extensive approach to dense, hot matter the particle and antiparticle oc¬ 
cupation numbers, rq and hj, can be obtained from the Jayne’s extremalization 
of the entropic measure 

S = i n i + (1 - rii) ln(l - m)] + [rii nj , (1) 

l 

under the constraints imposed by the total number of particles, N , and the 
total energy of the system, E (e* is the energy of the i-th energy level) [IB], 

(m - rii) = N and ^ (n* + n*) e* = E. (2) 

i i 

As a result we get the Fermi-Dirac distributions, 


Tli = 


exp (xi) + 1 


rii 


exp (xi) + 1 ’ 


which depend on the dimensionless quantities 


( 3 ) 


x = f3{e — fi) and x = f3(e + fi). (4) 

where (3 = 1/T, e = yfp 2 + m 2 , m is fermion mass and fi its chemical potential. 
For e = 0, the distributions n(x) and n(x) = n(x) satisfy following relation: 

n(x) + n(x) = 1. (5) 

The chemical potential (separation energy) for anti-fermions is negative — jjl 
(because the N increases after removing the anti-fermion), therefore this relation 
allows holes among the negative energy states to be interpreted as antiparticles 
with positive energy. As we shall see below, this is usually not the case for 
non-extensive distributions [IT]. The similar relation for particle distributions 

n(x) + n{—x) = 1 (6) 

exhibits, particle-hole symmetry around the Fermi surface where x = 0. 
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When formulating a non-extensive version of the model of dense matter one 
must do it in a thermodynamically consistent way. This means that one has 
to preserve the standard thermodynamical relationships among the thermody¬ 
namical variables such as, for example, entropy, energy and temperature, 

dS\ _ 1 
dEJv,N~r 



The best solution is to use once more Jayne’s maximum entropy prescription 
and extremalize the appropriate non-extensive entropic measures with specif¬ 
ically chosen constraints because, as was shown in [18], any thermostatistical 
formalism constructed by this method complies with the thermodynamical re¬ 
lationships [ |15j . A thermodynamically consistent formulation is then obtained 
by the appropriate identification of relevant constraints with extensive ther¬ 
modynamical quantities (like number of particles N or energy E ) and of the 
corresponding Lagrange multipliers with appropriate intensive thermodynam¬ 
ical quantities (like temperature T and chemical potential p.). In order to de¬ 
scribe non-extensive many-body systems characterized by q ^ 1, one has to 
introduce the non-extensive measure of the entropy expressed, in analogy with 
Eq. (CD), in terms of single particle distributions for fermions and anti-fermions, 
with accordingly modified constraints (ED - In this way we obtain a q-generalized 
quantum distributions from a given form of the non-extensive entropy. 


There are different choices of non-extensive entropy which lead to different 
thermodynamical solutions pdJlll5ll20j . We shall discuss three possibilities, two 
of which are currently exploited |9Fl5] . In the first choice, one uses a straight¬ 
forward (/-generalization of the entropic measure used in Eq.(Jl]), which now has 
the following form: 


S q = J2 s qi where for 0<n qi (n qi )<l 

i 

(\ 1 - n 9 qi-( 1 -n q i) q r v - , 

(•) s qi = - q-1 - K riqi^Tlqi}. (8) 

This generalization consists in replacing in Eq. (JTD ln(n) by In q (n q ) = (1 — 
n l ~ q )/{l — q) and in using as the corresponding effective occupation numbers 
of particles and antiparticles, n qi and n qi , instead of n % and n t . This reflects the 
fact that non-extensive properties can be formally understood as a presence of 
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some effective interaction between the constituents of the system, the strength of 
which is proportional to \q— 1| (in the non-extensive environment, parameterized 
by q, constituents are no longer free particles) [21]. Therefore constraints (ED 
for the total energy E of the system and th e co nservation of the total fermion 
number N will now have the following form 


irk 


E (n 9 qi - n qt) = N and E ( n h + n li) e* = E. (9) 


The extremalization of entropy given by Eq. (JSD , performed under constraints 
(P), with respect to two independent variables, n and h, gives the following set 
of equations with Lagrange multipliers a and f3: 

c) s c) s 

-^-=(a+(3e)qn i q q ^ 1) ] = (-a+/3e)qn < £~ 1) for 0<n q (n q )<l, (10) 

Solving these one obtains the following single particle distributions for fermions 
and anti-fermions: 


= with e «M = [i + (?-i)W, 

where x and x are the same as in Eq. (JTJ) with chemical potential fi = —a/(3. 
However, these solutions must be supplemented by conditions assuring that 
G q {x) in Eq. ( TTlT) is properly defined, i.e. that 

1 + (q — l)x > 0 

for all q and x{x). We get therefore four constraints for particle momenta p and 
antiparticle momenta p with q < 1 and q > 1: 



p 2 + m 2 < 


T 


q -1 

T 

p 2 + m 2 > - 1 - fi 

q- 1 


+ fi, \jp 2 + m 2 < 


T 


p 2 + m 2 > 


q -1 
T 

q -i 


— fi for q < 1; (12) 

— fi for q > 1. (13) 


These lead to the so-called Tsallis’ cut-off prescription for single particle distri¬ 
butions: 


3rk 


This form assures the fulfillment of basic requirements of thermodynamical consistency 119115120 ] . 
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( 14 ) 


n q< i = 0 

for e > fi + 

n q> i = 1 

for e < n — 


P 

( 1 -?)’ 

P 

(?-l)’ 


n q <i = 0 for e > 
no limitations for 


p 

^ (1-9) 

hg>l- 


(15) 


Notice differences in restrictions imposed on momenta p and p for q > 1 and 
q < 1 in Eqs. (TT21) and (TPl) . For q > 1 and sufficiently small momenta p, the 
limiting energy, eu m < —T/{q — 1) + p,, gives the usual prescription (TT5T) used in 
JSj. For q < 1 and for large momentap, the limiting energy, en m > — T/(q— l)+/i, 
cuts off large momenta, well above the Fermi level. However, limits given in ( TT2h 
for anti-fermion momenta practically prevent their propagation (as particle- 
antiparticle pair excitations) for small temperatures, see Eq.( lTD . Consequently, 
this choice of S q allow only for q > 1 in the propagation of virtual particle- 
antiparticle pairs. This means that for S q given by Eq. (JSD only the q > 1 case 
has physical meaning in relativistic Fermi gas model. Note that relations (El) 
and (ED are not now satisfied because e q (x ) * e q (—x) ^ 1. 


In order to allow for a relativistic calculation with q < 1 and including virtual 
particle-antiparticle pairs, we propose to use the entropy functional (cf. ED but 
with a different (dual) parameter q for antiparticles, 

q —>• q = 2 — q. 

This is the second choice to be discussed. In this case we have: 




v*/ 


Sqi 


1 - n ai ~ l 1 - U qi) q 


q -1 


+ for 0<n qi (n qi ) < 1 (16) 


As a result, after extremalisation, with the same conditions as before, one ob¬ 
tains the following occupation numbers: 


= ^ = MiKT’ with ^ ) = 11 + ( 9- 1)li ‘ 1 - (17) 

Note, that now the corresponding effective occupation numbers are n q q for par¬ 
ticles (as before) and n q q for antiparticles. The non-extensive analog of relation 
(ED is now satisfied because dual q = 2 — q in the anti-fermion distribution in 
Eq. flTTft means that e q (—x) = l/e q (x). Therefore 
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n q (x) + n q (-x) = 1. 


(18) 


To summarize, the q > 1 relativistic dynamics will be realized in the first 
choice of S q , given by Eq.flED, whereas q < 1 dynamics will be realized in the 
second choice of S q , Eq.( TT6j) . However, only in the later case is the property 
H81 satished. 

The third choice of non-extensive entropy consists in removing all previous 
additional cut-offs by using a new definition of entropy, S q (cf. ED- It has been 
proposed already in [T5] but without including anti-fermions. In our case this 
choice means that: 

(*) we assume q > 1 for fermions above the Fermi sea, 

(*) we assume q < 1 for fermions in the Fermi sea, 

(*) we assume q > 1 for anti-fernrions. 

The resulting s q i has the following form for q > 1 (cf. ED: 


(\ 


• 


• 

Sqi = “ 




L qi V - 1 ,L qi) 

< 7-1 


T \jX q i y Tl q i\ 


for 


1 n qi ^1 n q tf + _^} fo r 


Extremalizing S q with the same conditions as given by 


0 < n q i < 1/2 

(19) 

1/2 < n qi < 1 

Eq.(ED one obtains: 


l 

l+e 9 (/3(e-/x)) 

1 

1 +eq(P{e-n)) 

1 

1 + e q (/3(e + fi)) 


for 0 < n q < 1/2, 
for 1/2 < n q < 1, 

for 0 < n q < 1. 



( 20 ) 


In this case the corresponding effective occupation numbers are n q q for particles 
in the Fermi sea, n q q for particles above the Fermi level and n q for antiparticles. 

Because the q parameter is now different (dual) for particles below the Fermi 
surface and particles above the Fermi surface, the following non-extensive ver¬ 
sion of relation (ED (obtained similarly to (TT81)), which connects particle occu¬ 
pations around the Fermi surface is satished: 
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n q (x) + riq(-x) = 1 . 


( 21 ) 


Note that in a non-extensive Fermi gas one deals with effective occupations, 


n 


eff 


= K) 9 , 


as was mentioned before. Therefore, relations 


and 


for 


the single particle distributi ons n 


distributions n q ^ and n q ^. 


~k ~k ~k 


and n n will not be satisfied for the effective 


To summarize, in this work we derive the fermion n q (p) and anti-fermion n q (p) 
distributions in the Fermi gas model for 3 different choices of non-extensive en¬ 
tropy, which correspond to different choices of non-extensive parameterizations 
II5j . These choices result in different descriptions of dense NM [S1191I11U1311TT] . In 
the first two choices, given in Eqs.(jS]) an d (EES), the non-extensivity parameter 
q is constant for all values of n q (n q ). Limits imposed by Eq. fTTZl) on the anti¬ 


fermion distribution allows its propagation only for q > 1. Therefore, in the 
first choice the particle and anti-particle distributions can be obtained for the 
same q > 1 [5] but without particle-antiparticle symmetry (1 181). This symmetry 
is restored in the second choice, where we choose an anti-fermion distribution 
with dual q = 2 — q for anti-fermions for all values of n q {n q ). The third choice, 


Eq. 


introduces a jump in the q parameter at x = 0 for particles which 


results in discontinuous behavior in the single particle entropy, single particle 
energy and effective occupation, n e ^, which in the non-extensive case has the 
power q [131 • 


The first choice, Eq.flEj), realized for q > 1, with a cut off (TT5h for small momenta 
p, introduces in general positive correlations which increase the particle occu¬ 
pation numbers at the bottom of the Fermi sea. The second choice, Eq. (ll6l) . 
with cut-off (1151) for large momenta p in the tail of the Fermi distribution, in¬ 
troduces correlations which eliminate large momenta well above the Fermi sea. 
However, this scenario is not supported by recent experiments which exhibit 
strong short range correlations with a large tail in the distribution above the 
Fermi sea f22]- The third choice, Eq. (TH)fh involves a change of the parameter 
q in the particle distribution at the Fermi surface. This change of parameter q 
from (q = 2 — q < 1) to (q > 1) means that we expect different, non-extensive 
effects below and above the Fermi surface. However, the jump in the q value at 
the Fermi surface produces an energy gap near the Fermi surface. Whereas such 

* * ^However, the departures from unity are small because q ^ q = 2 — q in the power of the effective 
distribution riq . 
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a gap is well known in nuclear physics and is produced by strong pairing corre¬ 
lations in the nuclear ground state below the Fermi surface, it disappears when 
the temperature increases [23]- Therefore, the above non-extensive description, 
which assumes a change of the q parameter at the Fermi surface, acts in the 
opposite direction and will not describe this phenomenon. On the other hand, 
if we modify the change in q in such a way that it goes smoothly from q in the 
Fermi sea to q above the Fermi sea, within the finite range of momenta near 
the Fermi surface, then the occupation and energy gaps can vanish. In this way 
we will get the restoration of standard Boltzman-Gibs statistics, with q = 1, at 
the Fermi surface. 


Such jumps in q value are absent for antiparticles. One has to stress that their 
role is vital in NJL models. In the WM type of approach they are usually 
neglected because in these models the effective degrees of freedom are projected 
onto the positive energy states (nucleons). However, the negative energy states 
(anti-nucleons) cannot be neglected in NM [511241125] , In the relativistic mean 
held model there are possible corrections coming from the vacuum polarization 
which involve contributions from nucleon-antinucleon pairs [51155] . Therefore, in 
the non-extensive version of these models (like in gWM) [7118IfTT] . one has to 
decide whether higher order NN loops should be included and for q > 1 (first 
choice) or (q < 1) (second choice). The nuclear EOS described in quark and 
anti-quark degrees of freedom is discussed in the g-NJL model [4TJJ approach. 
Here the propagation of quark-antiquark pairs, described by q > 1 statistics, 
is essential for the dynamics of massless quarks which change their phase at 
the critical point. The non-extensive corrections spread the critical point over 
a wider area of temperature and density [TT] • 


A description of the EOS in terms of nuclear degrees of freedom is well for¬ 
mulated in WM [3]. There are two groups which generalize this mean field ap¬ 
proach. The constant non-extensive parameter q > 1 is applied in [TO]. The oth¬ 
ers [811111114] apply different distributions which cannot be obtained directly from 
extremalization but they satisfy basic thermodynamical relations [H]. Their n q 
distributions are formally identical to those obtained from our third choice (and 
also to those obtained in [501) but later they use 18111111141 a constant power q 
in n e /f = n q q . Only in [I5J the q power changes in the Fermi surface, like in our 
results for the third choice given in Eq.(TT9fb Therefore, only the results of [13] 
are fully consistent. 
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The following relation for the effective distribution n q ^ shows its behavior 
around the Fermi surface - where n q (x = 0) = 1/2: 


n q (x) + n q (—x) < 1 

for 

q > 1 

1st choice Eq.([8D, 


n q (x) +n q (—x) > 1 

for 

q < 1 

2nd choice Eq.(lTF|). 

(22) 

n q (x) + n q (—x) ~ 1 

for 

q « 1 

3rd choice Eq.(lT91). 



Here we assume for the third choice ( ITT)}) that the jump in parameter q is 
smeared for particle momenta around x = 0 with q = 1 at the Fermi surface. 
These properties should now be compared with realistic distribution of large 
nucleon momenta induced by short range correlations observed recently in [22] • 
Such a comparison should allow for a better understanding of the role of the 
parameter q. 

Concluding, different q parameterizations presented in this work predict dif¬ 
ferent non-extensive effects in different part of the particle spectra. The main 
difference was found at the Fermi surface between smooth non-extensive dy¬ 
namics with constant q , (with q > 1 for the first choice and q < 1 for the 
second one), against the behavior observed in the third choice (1T91). character¬ 
ized by a substantial change in the non-extensive parameter q across the q = 1 
value on the Fermi surface. Observations which will indicate the restoration of 
standard Boltzman-Gibs statistics {q = 1) at the Fermi surface would support 
the third ( TT9T) choice which is used [13j in the calculation of the EOS with nu¬ 
cleon and meson degrees of freedom. The relations ( 12011 which reflect the broken 
symmetries of extensive dynamics, should be accounted for in any non-extensive 
analysis in NM. It is interesting to note that changes in the non-extensive pa¬ 
rameter q (like assumed in (TT9h ) can produce energy gaps, which are in general 
connected with attractive correlations in Fermi systems. 
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